Ac hopping conduction at extreme disorder takes place on the percolating cluster 
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Simulations of the random barrier model show that ac currents at extreme disorder are carried 
almost entirely by the percolating cluster slightly above threshold; thus contributions from isolated 
low-activation-energy clusters are negligible. The effective medium approximation in conjunction 
with the Alexander-Orbach conjecture leads to an excellent analytical fit to the universal ac con- 
ductivity with no nontrivial fitting parameters. 



Recent advances relating to ion conduction in glasses 
and other disordered solids include the application of 
multidimensional NMR techniques Jj, the introduction 
of ac nonlinear spectroscopy 2], and elucidations of the 
high-frequency nearly constant loss 3- Moreover, it was 
found that the old idea of ions moving by the vacancy 
mechanism may well be correct [4] , and simulations gave 
new insight into the mixed-alkali effect Q . Despite these 
and other significant advances, important questions re- 
main unanswered. For instance, it is still not understood 

\ what role is played by ion interactions for the conductiv- 
ity @, or why the random barrier model (REM) 0, Q 

\ represents ac conductivity data so well. The latter ques- 
tion is not answered below, but new simulations and ar- 
guments are presented that we believe lead to a full un- 
derstanding of the physics of the RBM in the extreme 

\ disorder limit (low temperature limit). 

Ac conductivity is often studied also for amorphous 

\ semiconductors, electronically or ionically conducting 
polymers, defective crystals of various kinds, polaronic 

\ conductors, etc 0, Q. It is a longstanding observation 
that all disordered solids have remarkably similar ac con- 

' ductivities Universal features include : At low fre- 
quencies the conductivity is constant. At higher frcquen- 

\ cies it follows an approximate power law with an expo- 
nent less than one that increases slightly with increasing 
frequency. When measured in a fixed frequency range, 
the exponent converges to one as temperature goes to 

\ zero. The ac conductivity is less temperature dependent 
than the dc conductivity and obeys time-temperature su- 

\ perposition (sometimes referred to as "scaling"). The 

' frequency marking onset of ac conduction, cj,„, has the 
same activation energy as the dc conductivity. 

These and other observed features are reproduced by 
the RBM characterized [1, [l^ by five assumptions: 1) 
All charge carrier interactions including self-exclusion 
are ignored; 2) Charge carrier motion takes place on a 
cubic lattice; 3) All lattice sites have same energy; 4) 
Only nearest-neighbor jumps are allowed; 5) Jump rates 
DC exp(— iJ/fc^r) have random activation energies with 
distribution p{E). In the RBM the ac conductivity cr(w) 
relative to (t(0) as a function of a suitably scaled fre- 
quency becomes independent of p{E) in the extreme dis- 
order limit, i.e., when the width of p(E) is much larger 



than ksT [8|. Despite lack of non-trivial free parameters 
the RBM universal ac conductivity gives a good fit to ex- 
periment [8]; more refined models yield results that are 
close to those of the RBM [Tl| . 

It is well-known that the percolation threshold deter- 
mines the dc conductivity activation energy [l^ . At low 
temperatures the particles preferably jump across the 
lowest barriers. The highest barriers on the percolation 
cluster are bottlenecks dominating the low-temperature 
dc conductivity. If Ec is the highest barrier on the perco- 
lating cluster, one has a{0) ~ exp{—Ec/kBT) as T ^ 
[l2]. In order to have a non-zero dc conductivity of the 
percolation cluster, barriers slightly above the percola- 
tion threshold must be included. This defines the "fat 
percolation cluster" [^; on length scales shorter than 
its correlation length the fat percolation cluster appears 
fractal, on longer length scales it appears homogeneous. 

Understanding the RBM universal ac conductivity in 
terms of percolation arguments is much more challeng- 
ing. Traditionally the problem was approached 
"from the high-frequency side" by proceeding as follows. 
For (angular) frequencies w LUm there is a character- 
istic activation energy E{uj) < Ec for motion on time 
scales ~ 1/w; when uj decreases towards ujm one has 
E{uj) Ec- Links with E < E{lo) form finite low- 
activation-energy clusters. The cluster size distribution 
is assumed to determine the ac conductivity. Some time 
ago we proposed what amounts to coming "from the low- 
frequency side," namely that all relevant motion takes 
place on some subset of the infinite percolation cluster 
[SJ. Numerical evidence for this conjecture is given below, 
where it is shown that contributions from low- activation 
energy clusters outside the fat percolating cluster are in- 
significant. Moreover, it is shown that by assuming that 
not just a subset, but in fact the entire percolation clus- 
ter is important, an excellent analytical approximation 
to the universal ac conductivity with no nontrivial fit- 
ting parameters may be derived. 

The simulations of the RBM reported below refer to 
the Box distribution of activation energies {p{E) = I/Eq 
for < E < Eq, zero otherwise); ac universality in the 
extreme disorder limit implies that this distribution gives 
representative results [1]. The lowest temperature sim- 
ulated is given by /3 = 320 where (3 is the inverse di- 
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FIG. 1: Results for the ac conductivity at /3 = EalksT = 320 
in rationalized units @]; the frequency marking onset of ac 
conduction, ujm, is of order 10~^*. Ten independent 96 x 96 x 
96 samples were simulated, (a) Real part of oiuS) with two 
cut-off's: k — 3.2 and k = 12.8, averaged over the ten samples. 
The dashed line is the prediction of Eq. (2) empirically scaled 
to the k = 12.8 data, (b) Relative deviation from k — 12.8 as 
a function of frequency plotted for each of the ten independent 
samples. 



mensionless temperature, /3 = E'o/fcsT. For f3 = 320 
the jump rates cover more than 130 orders of magni- 
tude, making simulations quite challenging. We used a 
method based on solving the Laplace transform of the 
master equation numerically [14]. Conductivity data for 
(3 = 320 give an excellent representation of the univer- 
sal master curve for the RBM over the frequency range 
studied here [1]. 

In previously reported simulations ,^ we applied an ac- 
tivation energy cut-off above the percolation threshold, 
Ecut/Eo = Ec + k/P, where Ec = 0.2488 is the percola- 
tion energy for the cubic lattice and k a numerical con- 
stant. Jump rates for links with activation energies larger 
than Ecut were set to zero in order to be able to simulate 
large samples. Figure 1(a) presents the real part of the 
ac conductivity a{uj) — cr'{Lo) + i(j"{Lo) for k = 3.2 and 
k = 12.8 respectively ai (3 — 320. There is little differ- 
ence between the two data sets. The dashed line gives the 
prediction of the diffusion cluster approximation (DCA) 
combined with the Alexander-Orbach conjecture as de- 
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FIG. 2: Contribution from isolated clusters for the real part of 
the ac conductivity, ajci^-^), relative to cr' is^) as a function of 
the scaled real part of the conductivity, (j'{oj) = o' (lo) / a{Q) 
(cut-off: k — 6.4). The two dashed lines mark the relative 
masses of isolated clusters. Their contribution, however, is 
much smaller than their relative mass, showing that the dom- 
inant part of the ac conduction takes place on the fat perco- 
lation cluster. 



tailed below (Eq. ([2])). Fig. 1(b) gives the relative errors 
involved for different k values, taking k ~ 12.8 as repre- 
senting the "correct" data. As expected, the errors are 
largest in the dc regime and decrease with increasing k. 
Choosing k — 6.4 gives an error of just 1-2%. 

We proceed to investigate the behavior with k = 6.4 
in more detail. Applying this cut-off, the links with non- 
zero jump rate fall into two sets, the "fat" percolating 
cluster and all remaining finite isolated clusters. The 
latter do not contribute to the dc conductivity. Accord- 
ing to the traditional approaches based on cluster statis- 
tics, however, they give a significant contribution to the 
ac conductivity as soon as cj > ujrn 0; [13] ■ This was 
never tested numerically. Figure 2 presents the contribu- 
tion from isolated clusters a[Q{uj) relative to the full ac 
conductivity as a function of the real part of the scaled 
conductivity a = a{uj) /a{Q). The dashed lines mark this 
relative mass of the isolated clusters for jS = 160 and 
P = 320, respectively. The quantity a[^{uj) / a' (uj) , how- 
ever, is much smaller than the relative mass of isolated 
clusters for the range of frequencies covered in the figure, 
i.e., up to 10 billion times ujm (compare Fig. 1). More- 
over, for /3 — * CX3 the relative mass of isolated clusters goes 
to one, whereas we find that (T[fj{LL!) / a' (uj) is independent 
of temperature and stays insignificant. In summary, the 
dominant part of the low-temperature universal ac con- 
ductivity comes from the fat percolation cluster [8] with 
little contributions from isolated clusters. 

We now turn to the issue of analytical approximations 
utilizing the effective medium approximation (EMA) 
[13, El, El]. If G = Poit) exp{-iujt)dt where Po(t) 
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is the probability for a particle to be at a site given it 
was there at t = for a homogeneous system with uni- 
form jump rate, the extreme disorder limit of the EMA 
self-consistency equation is Intr = KfiibjG where A is a 
numerical constant [lo| . This determines a frequency- 
dependent complex "effective" jump rate that is propor- 
tional to the frequency-dependent conductivity 0, 1^. 
Henceforth we switch to the rationalized unit system [8| 
where the EMA selfconsistent ac conductivity equals the 
complex effective jump rate. Because Poit) is a function 
of the effective jump rate times time, crt, the quantity 
ibjG is a function of iuj/a. In the frequency range rel- 
evant for the universal ac conductivity of the extreme 
disorder limit corresponding to times obeying at ^ 1, 



one has \iujG\ < 1 0, [l7| . If d is dimension, whenever 
d>2 iujG as a function of iu /a has a regular first order 
term 0, i^G = ai{iuj/a) + .... If A is absorbed into 
a dimensionless frequency by defining Cj = aih(3uj / a{Q), 
the EMA universality equation 
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for d > 2 is 
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This equation gives a qualitatively correct, but numeri- 
cally inaccurate fit to simulations 8]. 

In our previous works it was proposed that some un- 
specified subset of the percolating cluster with fractal 
dimension d/ ("the diffusion cluster") is responsible for 
the ac conduction Q. If < 2 this led to the diffusion 
cluster approximation (DCA): Ina = {iuj / a)'^^ ^'^ [s*]. If 
the diffusion cluster is the so-called backbone, one ex- 
pects df = 1.7, if the diffusion cluster is the set of red 
bonds, one expects df = 1.1 [Sj]. Treating d/ as a fit- 
ting parameter led to d/ = 1.35 [8|], however, leaving the 
nature of the diffusion cluster as an open problem. 

What if not just a subset, but in fact the entire per- 
colating cluster contributes significantly to the universal 
ac conductivity? Random walks on a fractal structure 
are characterized by Po(t) oc {at)^'^"/^ [l^ where dn 
is the spectral dimension. For dn < 2 this leads to 
iujG oc (iuj/aY"^^ . In terms of a suitably scaled fre- 
quency the EMA thus implies the DCA expression with 
df = dff- According to the Alexander-Orbach conjec- 
ture [191] - known to be almost correct (see, e.g., 1J,120|) 
- one has dn = 4/3 for the infinite percolating cluster. 
If frequency is suitably scaled, this leads to the follow- 
ing approximation to the universal ac conductivity of the 
extreme disorder limit: 



Ina = I ^ 



2/3 



(2) 



As shown in Fig. 1(a) this expression provides an 
excellent fit to the universal ac conductivity of the ex- 
treme disorder limit [2l[. Equation ^ may be put to a 
more severe test, however, than just fitting the real part 
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FIG. 3: Testing Eqs. ((2| and ((3)1. Data represent averaging 
over 100 independent 96 x 96 x 96 samples (/3 = 320, cut-off: 
k — 6.4) and 100 independent 64 x 64 x 64 samples (/3 — 160, 
cut-off: k — 6.4). For both temperatures the frequency was 
empirically scaled such that a — 1 + uj in the low frequency 
limit, where a = a{uj)/a{0), uj = ujAeeo/a{0) Q|. (a) shows 
I In (t| as a function of \Cij/S-\ in a log- log plot, (b) shows the ap- 
parent exponent dln(| lna"|)/dln \i^/5-\ as a function of scaled 
frequency. A cross-over from fractal behavior (exponent 2/3) 
to homogeneous behavior (exponent one) is clearly visible, (c) 
shows the real and imaginary parts of the scaled conductivity 
compared to Eq.(|3l) (full lines). 



of (j{Cl'). Figure 3(a) tests one implication of Eq. 
Ino-j = Iw/o'l^/'^, by plotting |lnfT| as function of p/o"! 
in a log- log plot. A cross-over between two power-law 
regimes is seen, corresponding to a cross-over between 
Eqs. ^ and In Fig. 3(b) the apparent exponent 
dln(| ln(7|)/dln [a)/(T| is plotted as a function of scaled 
frequency. Similar results are found by plotting the ra- 
tio of the phases of the complex numbers In it and iuj/a 
(data not shown). 

The picture emerging from Figs. 3(a) and (b) is the fol- 
lowing: Equation ^ works well whenever (D ^ 1; here the 
fat percolation cluster appears fractal because over one 
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cycle the particles move less than the correlation length. 
At low frequencies there is a transition to the analytic be- 
havior predicted when the dimension is larger than two 
(Eq. ([1])); over one cycle the particles here move longer 
than the correlation length and consequently the fat per- 
colation cluster appears homogeneous. 

The entire frequency range is accurately described by 
the expression 



ln(7 = 



1-K2.66— 



'1/3 



that is plotted as the full lines in Figs. 3(a) and 3(c). The 
exponent — 1/3 was chosen to get agreement with Eq. ([2]) 
for \Cj/cf\ > 1. The difference between Eqs. ^ and ^ 
is significant only at such low frequencies that (7(a)) is of 
order unity (Fig. 3(a)). Equation ^ breaks down only 
for the imaginary part for uj < 1 where Eq. ([2]) predicts 
a" oc w^/^ instead of the observed a" (x uj. Numerical 
solutions of Eqs. ([2]) and ([3]) are provided Ref. [22 1. 

In our opinion, the RBM must now be regarded as 
solved in the extreme disorder limit in the sense that a 
good understanding of the model's physics is at hand, 
leading to an accurate description of the ac conductivity. 
A notable consequence of the above is that the EMA - 
generally believed to be inaccurate except at weak dis- 
order - works surprisingly well in the extreme disorder 
limit if the "geometrical" input G is taken to reflect the 
fractal geometry of the percolation cluster. It would be 
interesting to know whether similar results apply when 
the EMA is applied for the extreme disorder limit of other 
models. 

The centre for viscous liquid dynamics "Glass and 
Time" is sponsored by the Danish National Research 
Foundation (DNRF). 
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